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Abstract 

This work concerns the dynamical two-point spin correlation functions of the transverse Ising 
quantum chain at finite (non-zero) temperature, in the universal region near the quantum 
critical point. They are correlation functions of twist fields in the massive Majorana fermion 
quantum field theory. At finite temperature, these are known to satisfy a set of integrable 
partial differential equations, including the sinh-Gordon equation. We apply the classical inverse 
scattering method to study them, finding that the "initial scattering data" corresponding to 
the correlation functions are simply related to the one-particle finite-temperature form factors 
calculated recently by one of the authors. The set of linear integral equations (Gelfand-Levitan- 
Marchenko equations) associated to the inverse scattering problem then gives, in principle, the 
two-point functions at all space and time separations, and all temperatures. From them, we 
evaluate the large-time asymptotic expansion "near the light cone", in the region where the 
difference between the space and time separations is of the order of the correlation length. 
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1 Introduction 



The Ising model has a long and illustrious history, due to its simplicity of formulation, its 
integrability and its applicability to physical systems. We shall be considering the quantum 
model in one spatial plus one temporal dimension, a spin-half chain which describes the be- 
haviour of strongly anisotropic systems such as Dy(C2H 5 S04)39H 2 0, FeCl22H 2 and related 
compounds [8]. The Hamiltonian is written as below, with af' z representing the Pauli spin 
matrices at site i 



where g describes the strength of an externally applied magnetic field. For fixed J, the param- 
eter g controls the ground state of the system, with g < J favouring a state with spins aligned 
(spin up) or antialigned (spin down) along the z axis and g > J leading to a ground state with 
spins aligned along the direction of the magnetic field. At intermediate values, the magnetic 
field induces tunelling between the spin up and down states. There is a duality mapping from 
the g > J regime to the same model with g < J. The self dual value of the couplings, g c = J, 
is a quantum critical point. In the neighbourhood of this point, the low-energy, universal be- 
haviour of the model is described simply by the quantum field theory of free Majorana fermions 
with m oc g c — g (in conformal field theory language, it is the Ising minimal model deformed by 
the energy operator). We will study this region; it occurs around any quantum critical point 
in the Ising universality class, hence is more general than the quantum Ising model itself. 

The quantities that are of most interest and usually accessible experimentally in spin chains 
are the expectation values at temperature T of products of time-evolved spin matrices, such as 
Tr (e~ H '/ T al (t)cx|(0)) /Tr (e -/f/ / T ), on which we will concentrate in the present paper. These 
may be related to response functions to magnetic perturbations which can be measured by 
neutron scattering experiments, for example. The scaling limit of these quantities near the 
quantum critical point is obtained by taking oc J/\g — g c \, t oc \/\g — g c \ and T oc \g — g c \ 

as g — > g c from above (disordered regime) or below (ordered regime). 

In the scaling limit, these quantities are related in a non-trivial fashion to the expectation values 
of products of fields in the fermionic model, due to the non-local nature of the mapping to the 
field theory. They are correlation functions of twist fields related to the Z 2 symmetry of the 
fermions, as recalled in section 12.21 Despite considerable effort, no systematic method exists 
for determining the large-time asymptotics of these functions at finite temperature. The well- 
known mapping of the finite-temperature quantum field theory on infinite space to a vacuum 
theory quantized on the circle leads naturally to correlation functions with imaginary time. A 
lot is known about these correlation functions, but the analytic continuation to real time is 
plagued by singularities, so more work is needed in order to investigate the finite-temperature 
correlation functions in time-like regions. Recent approaches to the problem include a semiclas- 
sical method applicable to the regime of small temperature [2H [23], an approach, working in 
a similar regime, based on identifying the leading singularities of operator matrix elements [T], 
a virial expansion (in powers of soliton density) [22], and a low-temperature expansion using 
a regularisation of infinite- volume form factors and an appropriate re-summation scheme 
Here, we take a different approach based on solving a famous set of nonlinear partial differential 




(1) 
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equations (PDEs) satisfied by the thermal correlation functions [2H [HI d3] , a generalisation 
of the nonlinear ordinary differential equations (Painleve III) occurring at zero temperature 
[261 [23 [13]. The main element of these PDEs is the sinh-Gordon PDE which is integrable 
and possess a zero-curvature formulation, hence is amenable to solution by the classical inverse 
scattering method (see the books jT2], 0]). This is essentially a generalization of the Fourier 
transform to nonlinear PDEs. A linear scattering problem with spectral parameter is associated 
to the sinh-Gordon equation, such that the associated time-dependent scattering amplitudes 
(corresponding to the solutions of the scattering problem with incoming plane wave boundary 
condition at x — > oo or x — > — oo) are in correspondance with the solutions to the sinh-Gordon 
equation. The beauty of this approach is that the time evolution of the scattering amplitudes 
is simple. Hence, finding the scattering amplitudes at zero time (the initial scattering data) 
suffices to find them at all times. Then, there are known methods for mapping back these 
scattering amplitudes to solve the original problem. The difficulties are in calculating the ini- 
tital scattering data for a given initial sinh-Gordon field configuration, and in mapping the 
time-evolved scattering amplitudes back to solve the PDEs in a time-like region. 

In this work, we first calculate the initial scattering data associated to the thermal equal-time 
two-point correlation functions in the Majorana theory, by solving a Riemann-Hilbert problem. 
Our calculation is mainly based on the theory of thermal states and finite-temperature form 
factors, initiated in [9] in the context of massive integrable QFT (in the context of more general 
QFT, see [18j [2j [15]). It turns out that the scattering data are essentially given by one- 
particle finite-temperature form factors of twist fields, evaluated in [91 [10]. These scattering 
data are checked against numerical solutions to the scattering problem and found to be in 
good agreement. We then recall how these scattering data may be mapped to the two-point 
correlation functions at all times, via a set of coupled Volterra integral equations known as the 
Gelfand-Levitan-Marchenko (GLM) equations. The derivation is not entirely standard, because 
in our case, the solution of the sinh-Gordon equation has complex regions (hence, we really have 
a hybrid between the sinh-Gordon and sine-Gordon equations). This approach is compared with 
the known large-distance series expansion for the correlation functions at zero time difference, 
the form factor expansion in the quantisation on the circle, and shown to reproduce the series 
term by term for the first few terms. Finally, we solve the GLM equations analytically at 
leading order (i.e. up to "exponentially smaller" terms) in the region x + t^> \x — t\, m _1 , T _1 
(where x is the separation), with m/T ^ and m\x — t\ ^ 0, both in the time-like (t > x) 
and space-like (x > t) cases. This region is chosen because it illustrates how to "cross the 
light cone" in the GLM equations, and because it hasn't been investigated yet in the literature. 
We find that at leading order, the two-point function in the disordered regime is given by 
the analytic continuation of the one-particle contributions to the form factor expansion in the 
quantisation on the circle (that is, the one-particle contributions to the finite-temperature form 
factor expansion), up to a factor of an exponential in x and t that is beyond the reach of 
the PDE method. The advantage of this approach is that temperature and mass are generic 
and the set of approximations required to solve the GLM equations analytically are carefully 
controlled. We believe that the GLM equations we obtained can be used to obtain systematic 
expansions in any region of space-time, and should give precise numerical solutions. 

The paper is organised as follows: in section 2, we present the necessary background including 
a discussion of the theory of the Majorana fermion and twist fields at finite temperature, 
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the various Hilbert space constructions that we will need, as well as the partial differential 
equations satisfied by the two point correlation functions of twist fields in the finite-temperature 
Majorana theory and the associated linear scattering problem. In section 3, we determine the 
scattering data associated to the linear problem and compare to numerical solutions to the 
scattering equations. In section 4, we outline the derivation of the GLM equations which allow 
for the solution to the original problem to be recovered from the time-evolved scattering data. 
In section 5, we solve the GLM equations for equal-time case and for the region x + t ^> 
\x — t|, m -1 , T -1 . In section 6 we present the results for the correlation functions in the latter 
region. We conclude in section 7. The appendices contain details of the calculation of the 
scattering data and further discuss the derivation of the results in section 5. 



2 Background 

2.1 The Majorana fermion 

The theory of massive free Majorana fermions, with real anti-commuting fields ip, ip on two- 
dimensional space-time, is defined by the action (we will consider a positive masse m > 
throughout) 

A = % [ d 2 x (-i(>(d x + dt)il> + ${d x - <W " ™H>) ■ (2) 



Note that this action implies that the leading term of the operator product expansion (OPE) 
of fermion fields IS, as X — > ±t from the space-like region, 

,, (x , mm ~^— y fotWo.O)--^^. (3) 

We introduce here three quantisation schemes, or rather three Hilbert spaces on which the 
model can be defined: the Hilbert space of quantisation on the line, that of quantisation on the 
circle, and the finite-temperature Hilbert space discussed in [9j[l0]. The latter one is the most 
important for our purposes. 



2.1.1 Quantisation on the line 

Let's take the theory (J2J) on infinite space x G (—00,00), and quantise on equal-time slices 
t = constant. Then the Hilbert space, Ti, is a subspace of the space of field configurations on 
the line, and is constrained to form a module for the canonical equal-time anti-commutation 
relations 

{if){x,t),4){x',t)} = 5{x-x') , {4){x,t),i){x',t)} = 5(x-x') , {if>(x,t),$tf,t)} = . (4) 
The fermion operators evolve on successive equal-time lines through the equations of motion 

dip(x,t) = —4>(x,t) 

— Tfl 

diP(x,t) = -i/>(x,t), (5) 
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with d = \{d x + d t ) and d = ^(d x — d t ), and due to the normalisation chosen, they are hermitian, 
= ij; , ffi — -0. A solution is given by: 

^(x, t) = y d# e 0/2 (A(6) e ip ° x - iEet + A\0) e -*w*+«%«) 

0(x,r) = J A6e^ 6/2 {A(9)e ipeX - lEet - A\6)e~ ipeX+iEBt ) , (6) 

where the mode operators and their Hermitian conjugate v4^(#) satisfy the canonical anti- 
commutation relations 

{A\9),A{9')} = 5{9-9') , (7) 
(other anti-commutators vanishing) and where 

pe = m sinh 9 , Eg = m cosh 9 . 

Here 9 is the rapidity, and throughout this work we will also use the variable A = exp(0). The 
Hilbert space 7i is simply the Fock space over the algebra (JTj) with vacuum vector |vac) defined 
by .4(0) | vac) = 0. We will use the notation 

\9 1 ,...,9 k ) = A\9 1 )---A^9 k )\v ac ) . (8) 

A basis is formed by taking, for instance, 9\ > ■ ■ ■ > 9k- The Hamiltonian is 

/oo 
&9 m cosh 9 A\9)A{9) , (9) 
-oo 

and has the property of being bounded from below on 7i and of generating time translations: 

d d 
i[H, i/;(x, t)} = —iP(x, t) , i[H, $(x, t)\ = — V(x, t) . (10) 

2.1.2 Quantisation on the circle 

Now let's take the theory (j2J) on the circle of circumference j3, x G [0, 0\, (x — 0) = (x — (3). 
Quantising again on equal-time slices t = constant, the Hilbert space is a subspace of field 
configurations on the circle. Thanks to the Z 2 symmetry of the action (j2J) under change of 
signs of the fields, it is possible to define the fermion fields on a double covering of the circle 
and there are two sectors available for quantisation: periodic or Ramond (R), and anti-periodic 
or Neveu-Schwarz (NS) conditions around the circle. The Hilbert space is again constrained 
to form a module for the canonical equal-time anti-commutation relations (J4]), and fermion 
operators evolving on successive circles obey (0) and are hermitian. We will denote these 
operators by ijjp, ipp. With S = NS and S = R for the NS and R sectors respectively, and 
Zns = Z + 1/2 and Zr = Z, the solution is given by 

Mx,t) = -L V* _JilL= (a n e ipnX - iEnt + al l e- iPnX+iEnt ) 



ij p — 0n/2 

Mx,t) = --= (a n ,e^-^-ate-^ + ^) , (11] 

V2/9 n ^„ Vcosh0 n 
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where the discrete mode operators a n and their Hermitian conjugate a) n satisfy the canonical 
anti-commutation relations 

{a+, cv} = 5„y , (12) 
(other anti-commutators vanishing) and where 

p n = msinhe n = ^^ , (13) 
E n = mcosh9 n . 

Here we define the quantised rapidities 

, /27m \ . 
9 n = arcsinh — — . (14) 
\mp J 

The Hilbert space TC^ is simply the Fock space over the algebra (Tl2l) with vacuum vector |vac)f 
defined by a n |vac)^ = 0. We will use the notation 

|m,...,n fc )| = aj ll -"4jvac>| . (15) 

A basis is formed by taking, for instance, U\ > ■ ■ ■ > n^. The Hamiltonian is given by 

Hp = £s + mcosh^n a" n a n (16) 



and generates time translations: 

i[H s ^ fo(x, t)] = *) . i[Hp, M*> *)] = fcM^ t) . (17) 



In ([TBI) , we have included the vacuum energies Eg- For both sectors they are infinite, but their 
difference is finite and given by 



AS = S R - S NS = / d6 coshfl In 



oo 



oo 



j _|_ g— m/3cosh0 
^ g— m/3coshS 



2.1.3 Finite-temperature Hilbert space 

When considering a quantum theory at finite temperature, it is convenient to work with the 
finite-temperature Hilbert space (or Liouville space) of thermo-field dynamics [HI [2], [15] . For 
the present model (jSJ) on infinite space x G (— oo, oo), this Hilbert space was studied at length 
in [91 [10] . The main idea is to consider (a subspace of) the space End (TC) of operators on Tt, 
essentially the linear span of all operators formed out of products of any number of operators 
a{6) and cr(8) at different rapidities, and to endow it with the inner product (parametrised by 
the temperature T) 

TV (e.- M ' T A<B) , v 
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where A and B are any two elements of End (7i). In fact, like in the case of the quantisation 
on the circle, thanks to the Z 2 symmetry of the action (J2J), there is another inner product that 
can be used to provide a Hilbert space structure: 

A ' B ~ Tt(e-^W) ' (20) 

where U is the unitary operator that implements the Z 2 symmetry, Uip(x,t)U = —ip(x,t), etc. 
(note that it is normalised to U 2 = 1). It will be convenient for what follows to have a symbol 
for the infinite valuaj 

T ( p -h/t\ 

where V is the (infinite) volume. With S = NS for the Hilbert space £ NS with inner product 
( TTTTj) and S = R for the Hilbert space £ R with inner product (|2"U|) . we will use the notation 

l n = |vac £ ) , = |0)_, = \9) + , g s(^„ S( o,\ = \°> >-.+> 



as well as Dirac's usual bra-ket notation for inner products, where 

s = m 



d(0)={ 1 + eTEe ! (22) 

5 = R 

(we introduce here the shorter notation g±(9) for later convenience). This factor is included 
here because it provides nicer analytic properties to the matrix elements of twist fields that we 
will need below. Given an ordering of the rapidities, for instance 9\ > ■ ■ ■ > 9^ all such states 
form an orthogonal basis, with inner product 



5 < 



Any operator on Ti can also be seen as an operator on C s by, for instance, its left action on 
operators on Tt. We will denote the left action of any operator A £ C s by A £ End . In 
particular, it is easy to see that 

i(£)|vac £ > 5 = £ W> 5 , it^lvac^) 5 = (24) 

and it was shown in 0, [TU] that, for instance, 

A(9)\9') s + = g s _(9)\9, 9') s , + 5(9 - 9>)\ vaCc ) s . (25) 



1 It can be formally evaluated by evaluating the traces explicitly using (zero-temperature) multi-particle 
states on the line, with a momentum discretisation Sp = 2irV~ 1 where V is the volume, corresponding to 
(5(6*1 - 02 = 0) h-> m^cosh(9i /(2tt) 
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2.2 Order and disorder fields and their form factors 

As discussed in the introduction, the quantities of interest in the spin chain model (CD) are corre- 
lation functions of the spin matrices themselves. We will be interested in finite-temperature two- 
point correlation functions Tr (e~ Hl / T af (t)crJ(O)) /Tr (e _ff///T ) in the scaling limit. These corre- 
spond to two-point correlation functions of twist fields in the Majorana theory [IT]; twist fields 
are the natural local fields associated to the Z 2 symmetry of the theory {4>,^} —> {—i/), —^}- 
This correspondance is a result of writing, in the mapping from the spin model to the Majo- 
rana theory, pairs of spin operators, say at sites % and j, as exponentials of sums of bilinears in 
fermionic variables along the segment of the chain joining sites % and j. There are two types of 
twist fields: the order field a, and the disorder field \x. Correlation functions of the former give 
the scaling limit of the spin-chain correlation functions in the ordered regime, whereas those of 
the latter give correlation functions in the disordered regime. 



2.2.1 Definition 

As operators on the Hilbert space 7i of quantisation on the line, there are two representations of 
the order twist fields, which we will denote a + and er„ . They may be defined through equal-time 
exchange relations with the Majorana fermions (see for example [16]): 

^{x)a + {x') = {-l) e{x -' x,) a + {x')^{x) (x ^ x') (26) 
V>(x)a_(x') = (-l) e ( x '- x) a-(x')i;{x) (x^x') (27) 

(and the same holds for ijj \—> ij}), where 0(x) is Heavyside's step function. The order fields 
are further defined by stating that they have a minimal scaling dimension, and that they 
possess non-zero vacuum expectation values (as operators, they have bosonic statistics). They 
correspond to the case g < J in the scaling limit of (pQ), for the spin operator we are interested 
in. 

In correlation functions on euclidean space, with insertion of Majorana fermions, this defines 
fields that create branch cuts on their right (er + ) and on their left (cr_), when correlation 
functions are seen as functions of the positions of Majorana fermions. The shape of the cut 
can be changed without any effect on the correlation functions, as long as it does not cross 
other local fields that are affected nontrivially by the Z 2 symmetry transformation. At zero 
temperature, the cut can also be rotated without much effect, so that cr+ and a_ are simply 
related: 

(vacl^^i) • • ■ip(xj)a + {Q)ip(x 1 ) ■ ■ •^>(a4)|vac) 

= (-l) k (vac\^(x 1 )---^(x j )a40)ij(x' 1 )---ip(x' k )\vac) (28) 

(and the same equality holds for ip i— > ip) - note that j + k must be even for the correlation 
function to be nonzero. But such a rotation is more complicated at finite temperature, where 
in calculating traces, one considers states with both even and odd particle numbers k (= j). 

The disorder fields fi± are defined through the leading coefficient of the OPE's ipo± (as x — > 0): 

V(ar)ff±(0) ~ - , = =r P±(0) ■ (29) 
2\jixx + z0+ 
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As fields in euclidean correlation functions, they still have cuts in right/left directions, 
but as operators on the Hilbert space, they satisfy different exchange relations, because of their 
fermionic statistics: 

4>(x)ii + (x') = {-l) e{x '- x) ^ + {x')^{x) (x^x') (30) 
4){x)^{x') = (-1) 6( "V-(3:>W {x^x'). (31) 

Their leading OPE coefficients are 

^(x)fi ± (0) ~ - 1 a ± (0) . (32) 

The OPEs lead to the following hermiticity properties: 

o± = o± , a4 = ±A*± • ( 33 ) 

Note that operators with branch cut on the right and on the left are simply related to each 
other through pre- or post-multiplication by the unitary Z^-symmetry implementing operator 
U introduced in (1201) : 

Ua± = a±U = a T , = —fi±U = — /i T . (34) 

Also, we will normalise the operators according to the standard CFT normalisation 

(x ± (x,0)cx ± (0,0) ~ \x\~* , /i±(x,0)/i ± (0,0) ~ ±\x\~* . (35) 
Relations (1341) imply that 

a + (x,0)a_(0,0) ~ , fi+(x,0)n-(Q,Q) ~ U\x\-* (36) 

(note that these normalisations are different from those used in [9J [TO], but more convenient 
here for clarity). 

On a given sector 7i^ s ' R of quantisation on the circle, operators associated to the fields a, \i are 
ill-defined. Rather, these operators are maps from one quantisation sector to another. There 
are two such maps for every field, H NS -> H R and H R -> H NS , but we will denote both simply 
by a and /j,, as it will be clear from the context which map they represent. Essentially, the cut 
associated to these twist fields, in this case, should be seen as extending in (imaginary) time, 
hence changing the quantisation sector at early or at late times. 



2.2.2 Form factors 

The main objects from which asymptotics of the correlation functions we are interested in can 
be evaluated are matrix elements of twist fields between the vacuum and excited states, or 
form factors. We recall here the known results concerning these objects (below, O is a for even 
number of particles and \x for odd number of particles; other form factors are zero). 
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In the Hilbert space on the line, 7i, form factors were evaluated a long time ago [3] by solving 
a Riemann- Hilbert problem on the rapidity variables (and then afterwards by a variety of 
methods), and are given by 

(^c\Q ± ( X ,t)\e u ...,e k ) ={±)k . m n\\^U^r^ ft tanh(^). (37) 



(vac|cr±(0, 0)|vac) \27r 



l<i<j<k 



In the Hilbert space on the circle, 7i^ s,R , the form factors were evaluated much more recently, 
first using lattice methods [6j [7] , then using the "doubling trick" directly in the QFT model 
© [H], finally using the connection to finite-temperature form factors in that model [9j [TO] . 
Those that will be of interest to us, with the NS vacuum, are given by 

N g(va C |Q(x,t)|n 1 ,...,n fc )g r f] A£it+j:U ( x ^ itE ) 

N S(vac|a(0,0)|vac> R " ' 6 W 

where the function h + {9) is defined below, (j4TJ 





i<*<i<fc 



Finally, on the finite-temperature Hilbert space £ ' R , they were evaluated in [9], [JO] by deriving 
a Riemann- Hilbert problem on the rapidity variables, similar to the one for ordinary form factors 
but with important differences, and by solving it. Those that will be of interest to us, this time 
on the R Hilbert space, are given by 



(vaC £ |0 ± Qr,t)|fll,--- A)*,,,,., M± _ r|] =A£x+-£U e^-itEe.) 

R (vac £ |a ± (0,0)|va Cz: )R 1 1 

k 

nw*) n 

j=l l<i<j<k 



2 



where 



I in 

e 4 

h±{9) = — = exp 

V 2,71 



^ 2ri S wh(0-8>) 11 {l-e~ E o'/T 



(40) 



As explained in [9], finite-temperature form factors specialise both to the zero-temperature 
ones in the limit T/m — > (recall that T parametrise the inner product used on £ > R ) when 
choosing all signes 6j = +, as is expected, and also to form factors on the circle, through 
the relations (with time taken imaginary, t = —ir, for the operator on the finite-temperature 
Hilbert space): 

N ^(vac|0(-r, -zx)|m,...,n fc ) R = f| J R (vac £ |(9 + (x, -ir)\9 ni + —,..., nu + —>;,...,+ 




(41) 

with (3 = 1/T. Hence, finite-temperature form factors offer a way of evaluating form factors 
on the circle. In a sense, the finite-temperature Hilbert space is the "analytic continuation" 
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of the Hilbert space of quantisation on the circle towards positive imaginary eigenvalues of 
the momentum operator p — > iE [9]. In this sense, the vacua |vac,c) ' R should be seen as 
"lying at x — > — oo" and NS ' R (vac,c| as "lying at x — > oo". For instance, the branch cut of 
+ (x,ir) towards the positive x direction on the right-hand side of ( l4"TT) changes the boundary 
conditions in that region from periodic (coming from the R finite-temperature Hilbert space) 
to antiperiodic, as expressed by the NS vacuum of quantisation on the circle on the left-hand 
side, and it is important that this cut be in the "direction" of the vacuum R (vac,c| and not of 
the excited states. 

Note that (with all operators at x — t — 0) 

N |(vac|a|vac) R = R (vac|a|vac)Jf = R (vac £ |a ± |vac £ ) R Z~5 = NS (vac £ |a ± |vac £ ) NS Z5 = s T , 

(42) 

where [231 12Q] 
s T = 

ms2i 2 e S A 2 exp 



and A is Glaisher's constant (and recall that Z is given by (I2T|) ). The appearance of the 
constant Z is a consequence of the normalisation (I36I) and of the finite-temperature form factor 
theory developed in [9j [loJE 



[m(3f 




d0id0 2 



sinh #i sinh 9 2 



(2ir) 2 sinh(m/3 cosh ^i) sinh(m/3 cosh 6*2) 



x In 



coth 



61 — 6r 



(43) 



2.3 Finite temperature correlation functions 

Whereas zero temperature field theory is often a good approximation in particle physics, con- 
densed matter applications often require the use of finite temperature field theory. This presents 
a major technical challenge, since finite temperature correlation functions contain contributions 
from all states in the Hilbert space. Physical observables at finite temperature are statistical 
averages over correlation functions, with the Gibbs density matrix e~ H ^ T for a temperature T. 
We therefore consider traces like 

Tr(e^---) 

In order to calculate such quantities, one may use two Hilbert spaces: that of quantisation on 
the circle, and the finite-temperature Hilbert space. 

The first arises form the well known "mapping to the cylinder". With the euclidean signa- 
ture (i.e. in imaginary time), finite-temperature correlation functions of quantum field theory 
models in one infinite dimension are equivalent to vacuum correlations of the associated model 

2 The idea is that it is the two-point functions R (vac£|C'+(a;, t)0- (0, 0)|vac,c) R or 
N (vac£|C'+(a;, i)O_(0, 0)|vac£) N , with these choices of branch cut directions ±, that one must take to 
obtain the finite-temperature form factor expansion, hence the clustering on the finite-temperature Hilbert 
space. 

10 



(44) 



quantised on the circle of diameter j3 = 1/T. That is, finite temperature correlation functions 
may be written as 



(0(x, -it) ■■■) t = Ns,s>{e™ s l 2 ■■■) s<vac|(V-T, -ix) ■ ■ ■ |vac)f , (45) 

where s is the spin of O and Ms,s' is some normalisation constant. The operator Op{— r, — ix) is 
the corresponding operator in quantisation on the circle, with space variable — r (parameterizing 
the circle of circumference (3) and Euclidean time variable x (on the line). The sectors S,S' 
depend on the operators O, . . . inserted. The trace ( 1441) with insertion of operators that are 
local with respect to the fermion fields naturally corresponds to the NS sector, S = S' = NS, 
due to the Kubo-Martin-Schwinger (KMS) identity 

(0(x, t) ■■■) T = (-1)*> (0(x, t-i/3)-- -) T , (46) 

where fo is 1 for an operator with fermionic statistics, and otherwise. In this case, A/ns.ns = 1- 
This is the sector with the lowest vacuum energy. However, with insertion of twist fields, if the 
branches are chosen so that at x — ► +oo there is a branch cut, then we will have S = R, and 
if they are chosen so that at x — > — oo there is a branch cut, then we will have S' = R, with 
a factor Z~2 in the number J\fs,s' f° r every branch cut. In particular, if the unitary operator 
U is inserted, other insertions being local with respect to the fermion fields, then S = S' = R, 
with A/r,r = Z~ l and 14/3 — 1. From the viewpoint of the model ([T]), correlation functions 
correspond to choosing the branch so that the branch cut only extends between the twist fields, 
so that S = S' = NS. 



The representation (1451) leads to a large- distance expansion of finite-temperature correlation 
functions, through insertion of the resolution of the identity on the Hilbert space 7~tp S . We 
write it here for twist fields O = a or O = fi in the quantisation on the circle in the NS sector: 

N |(vac|O(-r,-za;)O(0,0)|vac)^ = ^ £ (47) 

k=0 ni,...,rii.£Z 

x N ^(vac|O(0, 0)K . . . , n k )* *(n x , . . . , n k \O(0, 0)|vac)£ s . 

The second way of dealing with finite-temperature correlation functions is using the finite- 
temperature Hilbert space. From the basic definition of this Hilbert space, it is easy to see 
that 

(0(x,t)---) T = NS (vac £ |a(x,t)---|vac £ ) NS . (48) 

Although the NS sector naturally appears here, it is possible, when twist fields are inserted, 
to use the R sector as well. Since we are interested in configurations giving a branch cut that 
extends only between the two twist fields O — a or O — fi, we may use 

NS (vac £ |C' + (x,t)C' + (0,0)|vac £ ) NS = Z' 1 R (vac £ |C> + (x, t)C>_(0, 0)|vac £ ) R . (49) 

This can be deduced from the relations (|34|) . One can then use a resolution of the identity on 
C R in order to obtain a large-distance expansion of the finite-temperature correlation functions, 
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for x > 0: 



R (vac £ |O+(x,t)O_(0,0)|vac £ ) R = (50) 

d0i ■ • -d6 k e -^+T.U^pe r itE 9j ) 

x 



E E 



x R (vac £ |0+(O, O)|0i, . . . , e k )l_ tk £l> ..., R (0i, • • • , 0*|6-(O, 0)|vac £ ) R . 



It is important, in deriving this expansion, that the cuts of the twist fields extend towards the 
"direction" of the vacuum states (see the discussion after ()4ip ); this is why we need to use 
the Ramond Hilbert space. It is also important that x > 0, although the region x < can 
be obtained easily by complex conjugation. By deformation of contours, this expansion along 
with (j32J) can be seen to reproduce (pT7|) [U [TU] . 

Both expansions above are convergent whenever x 2 > t 2 (and x > 0), with t = Im(r) in (|47p . 
The expansion (I50p may offer the hope of a series converging in time-like regions by deformation 
of contours, but this process is plagued by singularities, except for the one-particle contributions. 
In the following sections, we will show that the combination of the finite-temperature Hilbert 
space approach with the integrable nonlinear PDEs satisfied by twist-field correlation functions 
allows one to investigate time-like regions. 



2.4 Correlation functions and integrable differential equation 

The scaling limit of finite-temperature two-point correlation functions of spin matrices in the 
quantum Ising chain gives, depending on the regime, two finite-temperature correlation func- 
tions in the Majorana theory: 

G(x, t) = (cr+(x, t)a + (0, 0))t (ordered regime) (51) 
G(x, t) — (fJ,+(x, t)fi + (0, 0))t (disordered rgime). (52) 

At zero temperature, these correlation functions satisfy a well known set of ordinary differential 
equations [261 [23 El [13] . These equations can be generalised to partial differential equations 
at finite temperature [211 CHI [13] . The resulting PDEs satisfied by the correlation functions G 
and G can be expressed in terms of tp and \i defined as 

G = s 2 T e x/2 cosh(</V2) 

G = s 2 T e x/2 sinh(^/2) . (53) 
In terms of these functions, the PDEs are 

2 

— Ttl 

ddp = — sinh(2^) (54) 
8 

2 

— 771 

<9<9x = -g-(l-coslM>)) (55) 
d 2 X + (dp) 2 = B 2 X + (B V ) 2 = , (56) 
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which include the sinh-Gordon equation for ip. 

The sinh-Gordon equation (154j) may be rewritten as the following linear problem with spectral 
parameter A = e 6 [T2l H] 

(d x - A,)* = - A t )V = , (57) 
where ^>(x,t) is a well-defined 2-component vector- valued function of x and t, and where 

%_( 2id t ip m(Xe~ v - \~ l e v ) \ 

x ~ 4 V rn{Xef ~ X^e-v) -2id t ip J 

i ( 2id x ip -m(Xe~ tp + X~ 1 e ip ) 

1 " 4 \ -m(Ae^ + X~ l e~^) -2id x ip 

Well-definiteness of ty(x,i) for any initial conditions (and for any spectral parameter) implies 
that d t A x — d x A t = [A t , A x ], which implies ( 1541) . 

In fact, under evolution in real time, the function <p is expected to acquire complex values 
in time-like regions of space-time. This is a simple consequence of the fact that quantum 
operators do not commute when separated by time-like distances in relativistic quantum field 
theory. Hence, we really have a "hybrid" between the sinh-Gordon equation and the sine- 
Gordon equation. From the viewpoint of the PDE above, the complex nature of ip for a real 
initial condition at t = can be seen as a consequence of the singularity at x = of this initial 
condition. This singularity comes from the fact that \im x _> G(x, 0)/G(x, 0) = 1, and stays on 
the light cone afer time evolution. In order to properly define the solution to the PDE, we 
need a prescription for going through it. The natural prescription, to which our results will 
apply, is Wick's rotation (Feynman's prescription) Im(£) = 0~. This corresponds to the usual 
definition of real-time correlation functions as analytic continuations from those in imaginary 
time t = —ir, tgR, where they stay real for all real r (that is, where operators are always at 
space-like distances). 

Hence at imaginary times t = —it (with dt<p = id T tp), the function ip is real. We will make use 
of two symmetries enjoyed by the linear equation (d x — A x (x; 6))ty(x; 9) = (recall that we use 
A = e e ), leading to two transforms in the space of solutions, and \& t— > \E fV , with: 

A* x (x; 6) = -A x (x) * (x; 9) = 9) 
A x (x; 9 + in) = a 3 A x (x, 9)a 3 =>• * v (x; 9) = a 3 ^(x; 9 + in) . 

Note that only the second is a symmetry in real time. 




3 Determining the scattering data 

The classical inverse scattering method [T21 H] is a generalisation of the Fourier transform to 
the case of nonlinear partial differential equations with a zero curvature formulation. The 
sinh-Gordon PDE above has such a formulation, the linear problem, and hence is amenable to 
solution by this method. Let us consider solutions to the spatial part of the linear problem with 
incoming plane wave asymptotic conditions from x = oo (positive frequency), or from x = —oo 
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(negative frequency). It can be seen that these exist thanks to the fact that the function ip 
vanishes exponentially at x — > ±00. The two corresponding solutions are known as the Jost 
solutions, which we will denote by and form a linearly independent set of solutions that 
are analytic in lm(9) £ [0, 7r], which corresponds to the upper half plane of A. Let us focus on 
one of them: 



The functions a{9) (called the Jost function) and b{9) are known as the scattering data. Given 
any ip as function of x (with appropriately vanishing large- \x\ behaviours), the linear problem 
above fixes the scattering data. The classical inverse scattering method gives the inverse map: 
given the scattering data, a function </? of x can be determined. This is done by solving set of 
integral equations known as the Gelfand-Levitan-Marchenko (GLM) equations. Through the 
time part of the linear problem (1571) . the scattering data acquire a simple time dependence: 



Their time evolution gives, by the inverse scattering, functions <p(x,t) that solve the sinh- 
Gordon equation (1541) . The Jost function a{6) does not have poles in Im(#) £ [0,7r], since it 
is proportional to the Wronskian of the Jost solutions, a(8) = — det(\& + , \P_)/2. However, b{6) 
may have poles in that strip. The presence of these poles is related to the way in which the 
function tp vanishes as \x\ ^00. 

The solution to (154)) that corresponds to the correlation function should be uniquely fixed by 
requiring that (j)(x,t — iT" 1 ) = (j)(x,t) and by requiring appropriate asymptotic conditions for 
x — > 00 and x 2 — > t 2 . These conditions are at the basis of the Riemann-Hilbert problem that 
fixes the finite-temperature form factors evaluated in [9l[T0]. In this section, we will determine 
the initial scattering data that corresponds to the correlation function by direcly using the finite- 
temperature form factors. Then, in section HI we will explain how the correlation functions are 
obtained from the GLM equations. 

3.1 The special solution 

In order to obtain the scattering data, a pair of independent solutions to the spatial part of 
the linear problem are required (or at least both their x — > 00 and x — > —00 asymptotics). In 
this sub-section, we show how to evaluate both asymptotics at large \x\ of one solution which 
is related to objects in the QFT model. 

The objects in the QFT model leading to a particular solution to the linear problem (157]) are 
related to the following traces (recall the notation (144)) ) : 




(60) 



a(9,t) = a(6), b(6, t) = e iEat b{6) . 



(61) 



(62) 



(63) 
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It can be shown that 

* = * Sym = e"*/ V> ? 2 " ^ ) (64) 

is a solution to ( 157|) . with A = e e . The calculation is included in appendix [A], and follows closely 
that of p]. 

In order to solve its associated connection problem (i.e. the problem of finding the asymptotic 
form at x — > — oo knowing that at x — ■> oo, or vice versa), we can evaluate the large- \x\ 
asymptotics using finite-temperature form factors. Since only the initial scattering data are 
necessary, let us consider the case t = 0. Using a resolution of the identity on the Hilbert space 
£ R , just as in equation (|50|) . we have, as x — * oo, 

F(x, 0; 9) = Z- 1 R (vac £ |a + (|, o) ti{0) ^ (-|, o) |vac £ ) R 

~ Z' 1 R (vac £ |a+ (~ o) |vac £ ) R R (vac £ |i t (^) /t- (-|,o) |vac £ ) R 

= 2T 1 R (vac £ |a + (|,0) |va CjC ) R <?_(0) R (%_ (-- o) |vac £ ) R 

= -i4e- A£a: -^p_(e)/i_(0) . (65) 

Here, we use the functions (1221) and (|40~I) . In the first step, we inserted U 2 = 1 in (• • -)t and 
used ( )34|) in order to have the Ramond Hilbert space £ R , in the second step we inserted the 
first term of the resolution of the identity, which gives the leading large-x asymptotic, in the 
third step we used equation ([23$ and in the last step we used the complex conjugate of the 
form factor formula fl39l) for k = 1 and for the minus sign (with h±(8)* = h Zf (9)), along with 
the anti-hermitian property fl33l) for fi_ . It can be checked that inserting the resolution of the 
identity between the operators A* and using (125]) . gives the same result. 

The asymptotic form of F follows similarly, for x — > oo: 

F(x, 0; 9) = iT 1 R (vac £ | / i + (f , o) A^) a_ (-- o) |vac £ ) R 

~ 2T 1 R (vac £ |/} + (|,0) At(0)|vac £ ) R R (vac £ |a_ (-- o) |vac £ ) R 

= 4e- A ^ + ^^ + W^W- (66) 

For the asymptotics as x — > — oo, we take complex conjugates in order to have operators ordered 
form left to right with decreasing values of their positions (so that the finite-temperature form 
factor expansion can be used), 

F(x,0;9r^L + (~,0)A(9)cT + (^0)) . (67) 



2W v ' T V2' y/r 
Now this expression can be evaluated cLS X ^ OO, clS above 

ipg 



F{x, 0; 0)* ~ s 2 T e A£x+ ^- g_(6)h-(d) (68) 
which gives 

F(x,0;9) ~ 4e A ^-^ g_(0)h+(8) . (69) 
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Similarly, for x — > — oo, 

F(x, 0; 9) ~ is 2 T e A£x+1J ^ g + (6)h^ (6) . (70) 



In summary, putting back the (imaginary) time dependence, the asymptotic forms as \x\ — > oo 
of the solution coming from QFT is 

^sym = e 2 2 g + h + I \- e 2 + 2 ( ^ \ (71) 

X^-OO . lP0±_Egr f I \ _ipgx Egr I \ \ 

ie 2 2 ( x I -ze 2 + 2 ^_/^ + ( _ 1 J , (72) 

where the real exponential factor e ±A£x was absorbed into the function \- Note that this 
solution is not a Jost solution, but is the most symmetric one: 

^sym ^sym j ^sym sym ■ \'°) 



3.2 The other solution and the scattering data 

In order to calculate the initial scattering data corresponding to the correlation functions, the 
asymptotics of a second, linearly independent solution to the linear problem is needed. In this 
sub-section, we find these asymptotics and hence derive the scattering data. Inspired by the 
form of the symmetric solution of the previous section, let us look for a complete set of solutions 
with asymptotics (again, at r = 0) 

*, "3" «*F 9+ A + (j)-/J e - 1 F,.k_(_ 1 1 ) (74) 
* =~ ,,^ 9+h _ ( \ ) (75) 

* 2 (l-„) e *F 9+A+ (})-(l- / J) e -^/L( _\) (76) 

x->-oc . _!Pfi£ , / 1 \ 

-ze 2 ^_^ + I _ 1 I ■ (77) 

The quantities a and /? are functions of the spectral parameter 9 to be determined, and the 
asymptotics as x — > oo of the two solutions are chosen to be in agreement with the sym- 
metric solution of the previous section. The Jost solutions may then be obtained by a linear 
combination of the equations above, and the scattering data obtained: 

i{(3 - \)h- i(3g- . . 

a = j-n w— , b = — — . (78) 

{p - a)h + ((J - a)g + 

Since TrA^. = 0, the Wronskian det(\l/i, ^2) is independent of x. Using this, we find one 
constraint: 

13 = 1 + a. (79) 
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The symmetry-transformed solutions, \P and \I/ V , also satisfy the linear problem. Hence, 

- x ^oo . ipg* ( 1 

x->— oo * i^af , / 1 \ , / 1 

a e 2 I 1 J -0 e 2 I _ x 

and by the fact that \Px an d ^2 form a complete set of solutions, we find, from looking at the 
x — > — 00 asymptotics, 

i|r 1 = _ia**i — i/3** 2 . (80) 
This gives, from the x — * 00 asymptotics, 

1 = a *a + /r(l-a) , = a*(3 + 0*(1 - 0) . (81) 

Also, we have 

$1 x =°° -i/3(9 + m)e^g+h+ [ \j + M^ + i^e'^g-h. ( _^ 

S^-OO „W£ / 1 

-e 2 I _ x 

That is, ^\ = —i^2, from which we have 

1 - a(0) = 0(0 + ztt) , 1 - 0(0) = a(0 + ztt) . (82) 

Equations ((72J) and ([HI]) show that a£M (for ^ el), and (J7HD with (PJ) gives 

a(0 + ot) = - a (9) . (83) 

A large-0 analysis (see the appendix [B]) leads to the condition 

a(9) ~ 1 as -> ±00 . (84) 

The final constraint comes from the fact that a(0) is analytic in the strip lm(0) £ [0, n] - recall 
from above that this follows from the analyticity of the solutions ty± in this strip and from the 
Wronskian identity a(X) = — det(*+, \&_)/2. Using the analytic structure of the function h-(9) 
in the strip lm(0) £ [— 7r,7r] (with 9 n defined in (fl4j) where = 1/T) 

h- : poles = n + — , n £ Z + - 

Z7T 

zeroes = ra + — , n £ Z , 

the conditions above form a Riemann Hilbert problem for the function a(9) which can be stated 
as follows: 

• a(9) £ R for £ R 
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• a{9 + m) = -a{6) 

• a(9) ~ 1 as 9 — > ±00 

• a (9) has zeroes at 9 = y + 9 n for n G Z + | 

• «(#) is analytic for Im(#) G [0, 7r] except possibly for poles at 9 — y + # n for nfZ 
The unique solution is: 

a(g) = l-e~^ ^) = 1 _ e 2 -^/T = 2 ^W> (85) 

and the resulting analytic structure of the function a(9) in the strip lm(9) G [— tt, it] is 

a : poles 9 — 9 n ± — , n£Z 

zeroes 9 = 9 n ± — , n6Z+ - . 

These results have been shown to be in good agreement with numerical solutions to the linear 
differential equations, see figure [TJ 

Then, the scattering data flTHl) can be simplified further, using the identities h_(9)h + (9) = 
(27ra(6 l )) _1 (the functions h±{9) are related by complex conjugation, and also by h + (9 + in) = 
ih_(9)), to 

(recall the functions g±{9) and h±(9), fT221) and (T4U|) respectively). The form of the initial 
scattering data a(9), b(9) above constitute the main results of this section. It is clear that the 
Jost function a{9) does not have singularities in the strip Im(#) G [0, tt], as it should. It is also 
important to note that it does not have zeroes either in that strip. 

There is a useful check of these results, concerning the second Jost solution, defined by 

x^oo / 1 \ , _iP61 

W_ = ce 2 — de 2 



1 V -1 



X^-OO -iP6 x ( 1 

= e 2 



-1 

We see from the behaviour as x — > — 00 that it is simply proportional to ^2'- 



*- = ^r- ■ ( 87 ) 

9-K 



The coefficients are therefore 



i(l — a)q + . , i(l — 9 , 9 , , 

c = -i ^± = 2ig + = b(9 + m) , d = — — = -2rna 2 h 2 _ . (88) 

9- K 

The various Wronskian relations obtained from the pairs (\&+, ^-), (^+, ^-), (^+? , ^ r +) are 

d=-a , \a\ 2 + bc* = 1 , b* = -b . (89) 
Our expressions for a(0) and 6(0) may be readily shown to satisfy these requirements. 
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i/f+ upper 




Figure 1: Comparison of the scattering data to numerical solution. The numerical solution was 
obtained by first evaluating the function tp using its form factor expansion in the quantisation on 
the circle, (I47p . up to 6-particle contributions, then by solving numerically fl57|) from incoming 
plane wave behaviour at x — 10. (a) shows the real (red) and imaginary (blue) parts of the 
first (upper) component of the Jost solution \& + . This should be compared to (b), showing 
the difference from the scattering data prediction, which is supposed to be exact for x — > — oo. 
The simulation is for pg = 1 and m = T = 1. The results are in good agreement (to 5% for 
x = -10). 

4 The Gelfand-Levitan-Marchenko equations 

In the previous section, we showed how the scattering data, a(A) and b(X,t), are obtained, (1861 ) 
and (I6T1) (with A = e e ). In this section, we show how the solution to the sinh-Gordon PDE, 
(p(x,t), is recovered from the scattering data. We follow the derivation of [I], although we use 
slightly different arguments in order to deal with a generically complex tp(x,t). 

4.1 Fourier representations 

Throughout this sub-section we consider time-evolved quantities, but omit the explicit time 
dependences for ease of notation. The gauge transformed Jost solutions 

vf + = e~2 CT2i *+ and = e~2 CTz \I>_ 

have Fourier representation 

*l = e ipW2 ( \ \ + £ dy (Ul(x, y) + X^W^x, y^e 1 ™' 2 (90) 
C = e-^ 2 ( _\ ) + f dy (y 2 (x, y) + \~ l W 2 (x, y))e-^^ 2 , (91) 
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where Ui(x,y) and Wi(x,y) are two component vector kernels. The gauge transformed linear 
equation satisfied by \& = \]/ + and \& = is (o^ — A x )^> = with the connection 

A x (x; A) = , { f 2 ( Vl A -i?i) _ s J J = AAl + V*) + A . 

Integrals of the type J°° dyf(y)e zpxV are O (A _1 e* PA:c ) as A — > oo, as can be seen from using the 
equality 

4? 

< \-2gipA?y/2 _ e ip\y/2 _|_ g e ip\y/2 

m\ y 

and integrating by parts. Using this fact and that equality, we can equate terms of order A in 
the equation for ^ + to obtain 

A / i 

1 H Ai UAx,x) = -Aq(x" 



1 



This is equivalent to the following relation between derivatives of the field <p(x,t) and the 
kernels: 

(d t - d x ) ip(x) = 2 (Uiix, x) + - U[(x, x)") . (92) 

In the above equation, + and — indicate, respectively, the first and second component of the 
vector. Further, by equating terms of order A -1 and substituting equation (I9"2"j) . we obtain 

e Mx) = 1 + —Wx{x,x)- W^x.xY + — [U x {x,x)- UAx,x) + 

m m m 1 V / 

- ^ ( d *Ui& V) + + y)~) \ x=y ■ (93) 

Note that this differs from formula 13.21 in [I]. Hence, knowing the kernels (which are functions 
of time through the time dependence of the Jost solutions), we may deduce the function ip(x, t). 
The Gelfand-Levitan-Marchenko (GLM) equations express Ui(x,y) and Wi(x,y) in terms of the 
scattering data. 



4.2 The GLM equations 

The previous subsection showed how the kernels are related to the solution, ip(x,t), of the 
PDE. Here, we demonstrate how the kernels are obtained from the scattering data a(A) and 
b(X,t) which we have obtained in section I3~2"l The Jost solutions are related to their symmetry 
transformed versions by 

= -^rr^MV) + (94) 

a(X) a(\) 

* V _(A,t) = ~ 77T^ / +(A, t) + *M*_( A ,t) . (95) 
a{A) a{\) 

The gauge transformed version of the first of these equations is 

L -C(A,t) = *i(\,t) + b -t^lr + (x,t) , (96) 



a(A) _v ' y +v ' ' a(A) 
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into which we may substitute the Fourier representations of the gauge transformed Jost solu- 
tions, fl90l and ()9ip . It is convenient to define 



r(X,t) 



b(-X,t) 
a(X) 



(97) 



Following the convention used in [1], we will also make use of a definition of integrals on the 
real A line that avoids singularities: 



dA = lim 



dA 



dA 



(98) 



Multiplying (1961) by \^e ipxV ^ 2 for j E {0, —1} and integrating over A from — oo to oo with this 
prescription leads to the GLM equations: 



dz 



-—a z U 1 (x,y,t) = F (x + y,t){ x ) + j [Mv + z^U^x, z,t) + F_ x {y + z^W^x, z,t) 

(99) 

F_x (y + z, t)th(x, z,t) + F„ 2 (y + z, t)W x (x, z, t) 

(100) 

where use has been made of the following functions, for j e {0, —1, —2}: 

poo 



—a z W 1 (x,y,t) = F_ 1 (x + y,t)[ 



dz 



Fj(x,t) 



47T 
1 

Air 



d\X j e ip ^ /2 r(X,t) 
d\X j e ipxx/2 - iExt r{X,0) . 



(101) 



These GLM equations are defined only for y > x, and the value for y = x of their solution is 
defined by the limit y — > x + . 



5 Asymptotic solutions to the GLM integral equations 

5.1 The zero-time form factor expansion on the circle recovered 

In this section, we check that the machinery we have developed to solve the sinh-Gordon PDE 
for all time differences t, reproduces the known result for zero time difference. Recall that the 
functions Fj(x,t) in the GLM equations are integrals of the scattering data, which are simply 
related to the one-particle finite-temperature form factors. However, the complete t = form 
factor expansion contains contributions from all particle numbers and rapidities. It is a non- 
trivial validation of our method, therefore, that the solution to the GLM equations contains 
terms in the t = expansion with more than one particle. 
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First, let us re-write the functions Fj(x,0) as convergent sums instead of integrals. Since 
a(A) is analytic in the upper half plane, the poles in r(A, 0) contributing to the integral in 
equation fllOip for x > are a result of the dependence on b(— A, 0). In terms of 9 = log(A), the 
poles are at 9 — m/2 + 9 n for n G Z (recall the notation (fl4l)). The residues of the integrand 
at these poles are 

rexp(-mxcosh(fl n )/2) / f°° d9' 1 / l — e ~ mcosh ( e ')/ T \\ 

mcosh(0 n ) 6XP V 2tt cosh(0 n - 0') g V.l + e— -shW/T J J ■ t 102 ) 

Hence, the functions are given by the following sums 

2i e 9„-mxcosh(9„)/2 / /-co ^/ j / j _ g -m cosh(0')/T 

0) = S ™cosh(0 n ) 6XP V 2 J-oo 2^cosh(6 n -6>) bg (^i + e-^coshw/r 

j^ e -mxcosh(9 n )/2 / /-oo ^/ j / -y _ -m cosh(0')/T 

F_!(X, 0) = E 2 mcosh (^) ( 2 y_ oo 2^cOsh(e n -0') ^ (l + e-mcosh^/T 

^ y c -e„-7TCTCosh(e„,)/2 / f°° d0' 1 /l — e -mcosh(fl')/T > 

F_ 2 (x, °) = E mC osh(# n ) 6XP V 2 i^oo 2^cosh(fl n -#') bg V 1 + e— cosh (e')/T 



We can obtain a large- distance expansion for ip using the form factor expansion in the quanti- 
sation on the circle. With the functions 

F N (9 ni ,9 na ,...,9 ntr )= J] tanhf^-^V (103) 

n^^^„<~ at \ ' 



0<p<q<N 

and 



/ f°° d6' 1 , /l-exp(-m/?cosh0')\\ , r~n rr, , * 

the form factor expansion (1471) for the two-point correlation function of a operators is 



N ^(vac|cr(0, -ix)a(p, 0)|vac)^ s = (105) 

oo / N \ 

4e~ A£x E E [U\9^\ 2 )\F N (9 ni ,9 n2 ,...,9 nN )\^ 



n=o ni ,...,n N eZ \i=l 

N even 



The same expansion, but containing only contributions from terms with odd particle numbers 
TV, exists for the correlation function of /i operators. Let us gather terms into groups with fixed 
particle number and use a short-hand notation for them: 



G(x,0) =^(vac|a(0,-z.x)a(0,0)|vac)^ b = s^e~^ x (1 + 2 ff + A ff + . . .) (106) 

and 



G(x,0) = N ^(vac|/i(0,-zx)/i(0,0)|vac)^ s = s 2 T e~ A£x (l ff + 3 ff + . . .) . (107) 
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Then, ip may be expanded in terms of form factors via its definition in equation (1531) as 



<p(x, 0) = tanh 1 



G(x,0) 



G(x,0) 

This is naturally a large distance expansion. 



[ ff + ( -3- - '//•-// + 3 /.f ) " 



108) 



In the GLM equations, the functions Fj(x,0) contain terms which are exponentially small at 
large distances x, as is clear from the sum expressions above. The iterative solution to the GLM 
equations, which is an expansion in powers of the number of functions Fj(x, 0), is therefore also 
an expansion of the type of the form factor expansion. It is natural to conjecture that the 
number of particles in the form factor expansion is in correspondance with the number of 
powers of Fj(x,0) from the iterative solution of the GLM equations. Equation f[§3"j) may be 
written as 



<p(x, 0) = - ln(l + I glm + Zglm + 3 GL M + • • •) 

1&LM / % GLM l&LJvA f^GLM ^GLM^GLM ^GLM 



4 J \ 2 2 6 



(109) 



in which 2qlm, for example, contains all terms present on the right hand side of equation (l9~3l) 
with quadratic powers of the functions Fj(x,0). Our conjecture is that expansions (11081) 
and (I109p should agree term by term. Let us show this explicitly for the linear term. The 
associated contributions to equation fl93|) containing only single powers of the functions Fj(x, 0) 
are 

— Ai - — - 4d - — - 16 / - — - - — - \ 

Iglm = W x (x, x, 0)+ + -Wxix, x, 0)~ - — (d^x, y, 0)+ + d y U x {x, y, 0)" \ y=x , (110) 

m m mr V ' / 

in which all terms in U(x, y, 0) and W(x, y, 0) are evaluated to one power of the functions 
Fj(x, 0). This is equivalent, therefore, to neglecting the contribution from the integrals, since 
these are quadratic or higher in powers of Fj(x, 0). We obtain 

iT e -mcosh(e n )x / roo d Q, j / j _ & - m cosh(0')/T \ \ 

^ 1(X ' X ' 0) cosh(g w ) eX H 2 |oo^ cosh(g w -gO lo ^ l + e— M^ jj 

Wt(x,x, 0)~ = -Wl(x,x, 0) + 

J" e -mcosh{0 n )(x+y)/2 / /•oc ^qi j /-^ _ e -mcosh(0')/T\ \ 

Z7l(x ' y ' ° )+ = ^ ^ ^ cosh(^) 6XP ( 2 2^ cosh(6> w -6K) l ° g { TT^^ J J 
th(x,y,0)- = -U 1 (x 1 y 1 0) + . (Ill) 

Then, only the Wi(x,x, 0) terms contribute to Iglm, since the contributions from t/i(x, y,0) + 
and Ui(x,y,0)~ cancel. We find 



^rj] g— m cosh(6 n )x / poo ^Q/ /j g— mcosh(S')/T 



cosh(# n ) 6XP V 2^cosh(# n -fl') ° g U + e" mcosh (e')/T 
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in agreement with Iglm = 4 Iff, with the identification 1/(3 = T. We have also shown 

that ^ 2 glm _ 1 glm \ = o by interating the integrals in the GLM equation once, and also that 

agreement is obtained for the cubic terms in the expansions. The three particle form factor, 
3//, contains the terms | tanh ^'"^ )| 2 . The fact that the GLM equations recreate these terms 
is a non-trivial validation of the method, since the building blocks of the equations, Fj(x,0) 
contain no such hyperbolic tangent functions themselves. 



5.2 Large-time asymptotics 

We now seek a solution that is valid in the region t + x ^> t — x, T _1 , m _1 , with t > x > 0. 
Our solution will be valid for m/T ^ and m(t — x) ^ 0. In that region, the residues of 
the integrand in equation ( 110 II) are no longer exponentially damped away from 9 = 0; a new 
approach is required. The main problem is that we need to cross the light cone x 2 = t 2 . The 
choice of the region above, "near" the light cone, is for two reasons: it is the simplest region 
where we can illustrate how to cross the light cone in the GLM equations, and this region has 
not been investigated in the literature yet. 

Notice that as A — > oo, we have r(A, 0) ~ 2%. Hence, from (UOip . there is a delta- function con- 
tribution at x = 2t in Fq(x, t). This is simple to extract, and we can write Fq as a distribution: 

2% 

F Q (x, t) = —5(x - 2t) + F P (x, t) . 
m 

The function Fq (x, t) is the principal part of F (x,t): it is equal to F (x,t) for \x — 2t\ > e 
and is zero (or any finite number) otherwise, where e is some positive number that needs to 
be taken to zero after integration over x. Taking the principal part is important, because the 
asymptotics of r(A, 0) also says that a pole will be present at x — 2t. Hence, the integration 
over x of an integrand that has F p (x, t) as a factor is the usual Cauchy principal value integral 
about that point. The delta-function contributes to the GLM equations in the time-like regime 
only, t > x. Let us consider the first GLM equation, (|99|) . The kernel U\ also acquires a 
delta-function piece, and we can also write it as a distribution: 

Ul(x,y) = 5(x + y-2t) ( ~ l ) +C(x,y,t), 

where that £ is a two-component vector that does not have delta-function contributions. It is 
expected to be itself a principal part of U\ with respect to various points, including x + y = 2t. 
The first GLM equation can then be written as (in the regime t > x): 

^a z £(x, y, t) = F p (x + y, t) (\ J + F Q p (y + 2t-x,t)(~*j 

+ J [F p (y + z,t)((x,z,t) + F_ 1 (y + z,t)Wi(x,z,t)jdz+^(x,2t-y,t) ( j \ 0(2* - x 

(112) 

where is the step function. This should be seen as an equation for a distribution £(x,y,t), 
and integrals as principal value integrals. 
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In appendix [UJ we give a derivation of the asymptotic solution to these integral equations, up 
to the integral over F_iW\. 

m f -1 



Z(x,y,t) = -Ff(x + y,t)l 1 j . (113) 

Here we will just check that this function gives equality in (11121) up to the integral over F_iWi, 
and up to terms 



V*M*+i0 2 /4-io+ (( m t)-~) , t - x, t - y fixed. (114) 



It will then be checked that the integral over F_iWi gives also this type of term. The derivation 
given in appendix [C] starts from a more general and natural ansatz and shows that the solution 
(11131) is the only one possible. We will not go into discussion of the uniqueness of the solution 
to the GLM equations; uniqueness proofs can be found, for instance, in [12] for the sine-Gordon 
model, but we do not know if uniqueness has been assessed in our case, where ip has a complex 
region. Naturally, the solution for </? with appropriate initial condition at t = is unique (with 
the extra condition that <p(x,t) = tp(x,t — i0 + ) in order to go through the singularity on the 
light-cone x 2 = t 2 ), and our solution to the GLM equations should give it. 

First, in the region x/2 + t 3> x/2 — t, the integrals in equation (1 10 II) are dominated by 
contributions from large 9 (i.e. as A — > ±oo). We can then expand r(6,0) and r{9 + in, 0) in 
powers of e~ e . The integrals for Fj(x, t) may be done analytically on each term of the expansion, 
and give an expansion in that region. We have 

oo 

r(9, 0) = 4ng + (9)h + (9) 2 = 2i ^ c^{T)e^ e e~^/ T (115) 

H,v=0 
oo 

r(9 + m,0) = -4ng^{9)h^{9) 2 = 2i ^ c^(T) e -^e- vEs/T , (116) 

ix,v=0 

with coefficients, c^ u (T),c^ u (T), coming from the expansion of the functions (I2"2"j) and (1401) in 
powers of e~ e . They are, in general, integral expressions with dependence only on temperature 
(and mass). Note, however, that c 00 = 1 and c^o = for all \i. With use of the identity 



oo 



u c e -au-b/u du = 2 ( ^ j 2 K_ 1 _ C (2V / ^) 



the above form of r(9,t) leads to the following expression for the functions Fj(x,t): 

Fjix, t) = ^jt cA T ) ( t ~ 2f 2 + Wr ) 2 K ^+j-^ m Viit + W2 + v/T){it - ix/2 + v/T)) 

fj,,i/=0 ^ ' ' ' 

i / m\ f —it + ix/2 + v/T\ 2 . n ; tt^t, ; tt^Tn 

- c ^ T ) { -it-J/2 + J/T j Ki+j-MVHt ~ ix/2 + v/T)(-it + ix/2 + u/T)) 

(117) 
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where K s (x) is the modified Bessel function of the second kind. Above, we take the principal 
branch of the square-root and in the cases where v = 0, we consider v = + . For j = 0, the 
above gives Fq (x, t) for x ^ 2t. Note that the term with \i — v — in the first series gives the 
expected pole 2t. 

In order to obtain terms up to (11141) . we disregard all terms with v ^ 0; in particular, the second 
series above can be neglected. Now we see that on the right-hand side of f II X 2 [) . the integration 
over Fq(u + z, t)£(x, z, t) has one or two poles about which the principal value has to be taken: 
one at z = 2t — y if y < 2t — x, and one at z = 2t — x coming from £(x, z, t). These principal 
values can be evaluated by closing the ^-contour around these points, and substracting the 
resulting half-residue. Closing the contour can only be made in the upper-half z-plane, because 
of the branch structure of the argument of the modified Bessel functions in the first series in 
the expansion for Fo above. Doing this exactly cancels the terms F^(y + 2t — x, t)(—i, i) T and 
— 2t — y, t)(l, l) T Q(2t — x — y) that were coming from the delta-functions. We are left with 
an integral of the form 

/ F (y + z,t)F (x + z,t)dz (118) 
J c x 

where the contour C x starts at x and finishes at oo by going through the upper-half plane. 
Everywhere on this contour, the leading asymptotic iT 1 _ /i (it) ~ \/^J, e ~ U °f ^ ne modified 
Bessel functions at large argument can be replaced, and the resulting integrals are conver- 
gent. These integrals are of the type (I114p . so can be neglected. More precisely, they are 
exp[— im\Jt 2 — x 2 — im^Jt 2 — (x + y) 2 /4]0(t a ) for some power a. For the case p* — 0, there is 
a logarithmic singularity as y — >• 2t — x, but no other singularities occur. Hence the neglected 
terms in £(x,y,t) can have at most logarithmic singularities, which are integrable, so under 
integration in (11121) they give again terms of the same type (JTHJ) • 

Hence, we have 

. YTl ( 1 \ 

Ux{x,y,t) = — F {x + y,t) ( 1 J + terms of the type (EI]). (119) 

This gives no contribution to tp(x,t) to the order of interest in equation (J93l) . since it is a 
function of x + y and proportional to (1, — 1) T , and the quadratic terms in equation ( 1931) can be 
neglected. On substitution into the GLM equation for W\, excluding the term corresponding 
to the integral over F^Wi, this result leads to: 

^* g Wi(x,y,t) = F_i(x + j/,t) ( M + J F_ 1 {y + z,t)F {x + z,t)dz. (120) 

Again, the integral part gives no contribution to the leading terms (since the delta-function 
piece is cancelled by closing the contour of the principal value integral in the upper half plane), 
so that we have 

j i \ 

Wi(x 7 y,t) = —F^x + ytt) ( _ x )+ terms of the type (HUD. (121) 

Recall that we have as yet ignored two of the integrals in the GLM equations: j F_xW\ and 
/ F_ 2 Wi. The results above vindicate this omission; such integrals lead to terms of the type 

dm. 
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Note finally that the space-like region x + 1 ^> x — t with x > t > can also be investigated 
with similar techniques, with the simplification that no delta-function contribution appears and 
that no singularities have to be dealt with. The results for U\ and W\ above are still valid in 
that region. 



6 Results 

All our results are valid for \t — x\m 7^ and m/T 7^ 0. The results of the previous section can 
be combined with (193]) to give 

tp(x, t) = -2iF_i(2x, t) + e ~ im ^^O ((mt)-°°) , t-x fixed. (122) 

That is, 

<p( x ,t) = -V cJT) I — ^] KMmVt 2 - x 2 ) + e~ imVW ^0 ((mt)" 00 ) , t - x fixed, 

(123) 

where c M (T) = c^o(T). Recall that the coefficients c^q(T) come from the expansion of the one- 
particle finite-temperature form factors in powers of e~ e . They may be calculated as integrals 
over rapidites, and are given by the generating function 

00 

The expansion f!123j) is a solution to the linearized version of the sinh-Gordon equation (1541) 
(that is, the Klein-Gordon equation). For convenience, we will use here the variables 

v = t — x , w = t + x . 

In fact, (I123p is an infinite linear combination of solutions w) = (v /w)% K^im^/vw) with 
spins [i — 0, 1, 2, . . . and temperature-dependent coefficients c M (T). 

In order to obtain the correlation functions, we need to evaluate x through (|55|) . (|56|) . These 
equations for \ leave room for three undetermined, generically temperature-dependent con- 
stants: 

X = -(A + Bx + Ct) + x (125) 

with lim^^oo x = 0. At T = 0, the constants A, B, C all vanish. In general, they cannot be 
determined using the integrable PDEs for correlation functions; an alternative method would 
be required. We expect that B + C > 0, which represents a decay of the thermal correlation 
functions at large time in the region we are looking at. 

The leading terms for x are obtained by keeping only the quadratic power of tp in (155]) . It is 
a non-trivial fact that any solution $^(1;, w) to the Klein-Gordon equation gives rise, through 



2i 



71 ^ 

w=0 



,2u)+l 



oo-i0+ / i 1 -mcosh(0)/T 

d6e^ e In ( i±f _ 



(124) 
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(j55p and (1561) . to a consistent large-?/; expansion for x of the form e 
growing integer powers of w~^. 

The large-w expansions of <p and x can be written 



-2imv 'vw 



times a series in 



<P 



X 



-2i 



urn 



vw 4 e 



2ixm 



vw 2 e 



-1vm*Jww 



1 + 



m^/vw m 2 vw 
fi(v) f 2 (v) 



From (11231) . we immediately have 



03 W 



m^^vw) 2 
/a(«) 



m\pvw m 2 vw m 3 (vw)2 



+ 0((mw)- 2 ) 
+ 0((mw)- 2 ) 



(126) 
(127) 



— h Cimu 
8 

9 3ici 



128 
75i 



mv + c 2 (mt> )' 



15ci 



15zC 2 



[mv) 2 + C3(mv) 



1024 128 

Then, solving (1551) to quadratic order of 99, or solving (1561) . we find consistently 

3z 

fi(v) = — + 2c 1 mv 

c / \ 33 ic~^ , o « , . q 

W) = + ^mv + {2c 2 + c^imv) 

/a W = ~ ~\§ mV ~ 4^ Cl + 10c Vy mv ) + 2(cic 2 + c 3 )(mw) 

The correlation functions (1531) are obtained from these series expansion: 



G 
G 



*-A-Bx-Ct , j_ 



2 



2 -A-Bx-CtP 



2 +e 



+ e~ 
0((mw) 



'0((mw)-°°) 



For G this gives: 



G 



-A-Bx-Ct 



1 + e - 2im v^ X 



7z + 8cimf 117 — 48icimt> — 32 (cf + 2c 2 )(?m>) s 



(128) 



(129) 



(130) 



(131) 



]irm 2 vw 327rm 3 (vw)2 



25Qirm i (vw) 



0((mw) 



It is interesting to see that the first two terms inside the parenthesis are temperature- independent. 
Taking T = (where A — B — C — 0), they can be seen to agree with the standard zero- 
temperature form factor expansion, which is valid also in the time-like region. The next terms 
acquire temperature corrections. Note that it was necessary to calculate and in order to 
obtain the last term, although it only involves the constants c\ and c 2 . 
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Note also that the same series expansion holds in the space-like regime x + 1 3> x — t, T -1 , m _1 
with x > t > 0. In this case, we need to take v with a phase e~ l7T . Since there is a path from 
v > to v < through the negative imaginary w-plane on which the series expansion is valid 
at every point, we expect that the constants A, B, C are the same in this space-like regime. 

From (11301) . the expansion (11261) can be used directly to give an expansion for G, valid for 
x + t ^> \x — i|,T _1 ,m _1 with x — t positive or negative (but, again, non-zero). In fact, the 
formula 

G = -is 2 T e- A - Bx - ct (F_i(2x,t) + e-^^^O ((mi)" 00 )) , t-x fixed 

shows that it is, up to the pre-factor e - A - Bx -ct ^ j ug |- ^ e finite-temperature form-factor ex- 
pansion up to one particle. The one-particle finite-temperature form factor expansion, initially 
an expansion at large distances, can indeed be continued to a large-time expansion without 
problems (in particular, to the region we are looking at here). Of course, the higher-particle 
contributions, under such continuation, would give contributions at the "one-particle" order as 
well, because of the poles in rapidity space (as briefly discussed in section I2T31 . Our results 
show that in our region of interest, these contributions only enter the exponential pre-factor, 
determining the constants A, B, C. In fact, the constant C is certainly expected to be modified 
in going from the large-distance expansion to the expansion in the region we are looking at. 
Indeed, in a large- distance expansion it is 0, but there is no reason to believe that it should be 
zero here. 



7 Conclusion 

Inspired by the link with finite-temperature correlation functions of spin operators in the trans- 
verse Ising model near its quantum critical point, we have investigated the dynamical two-point 
correlation functions of twist fields in the free Majorana theory at finite temperature and real 
time. These objects are hard to calculate by the well-known mapping to the vacuum theory on 
the circle, since this approach leads naturally to imaginary time and the analytic continuation 
to real time is generically plagued by singularites. Instead, we used the fact that the correlation 
functions satisfy a well-known set of integrable partial differential equations (PDEs), including 
the sinh-Gordon equation. These are amenable to solution by the classical inverse scattering 
method, which is a generalisation of the Fourier transform to tackle such nonlinear, integrable 
PDEs. The method relies on solving the associated scattering problem and mapping the results 
back to the solution of the PDEs. The main objects are the scattering data, related to scat- 
tering solutions with incoming plane wave boundary conditions at ±oo. The scattering data 
evolve in a simple way with time, and an important step is to determine the initial scattering 
data. 

We evaluated the initial scattering data that correspond to the correlation functions at finite 
temperature using the Liouville space of QFT (Hilbert space of finite-temperature states), and 
found that they are simply related to the finite-temperature form factors of the twist fields. 
Following the classical inverse scattering method, the solution to the PDEs are then obtained 
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from the solution to a pair of coupled Volterra linear integral equations (the GLM equations). 
These equations (I99lll00p . with the kernels (110 II) containing the initial scattering data, and the 
relation to the sinh-Gordon solution (1931) . are part of our main results. We showed that these 
integral equations reproduce the known results for finite-temperature correlation functions at 
zero time (the form factor expansion in the quantisation on the circle) term by term for the 
first few terms. From the GLM equations, we then evaluated the expansion of the correlation 
functions at large times, in the region where the space variable is kept of the order of the 
time variable. To our knowledge, this region, "near" the light-cone, had not been investigated 
before in the literature. We found that the one-particle contributions to the finite-temperature 
form factor expansion, which naturally gives an expansion at finite time and large distances, 
also gives the right expansion in that region, up to exponentially smaller terms, and up to a 
factor of an exponential in space and time which is beyond the reach of the present method. 
This is non-trivial, since higher-particle terms in the form factor expansion are expected to 
give contributions to the same order as that of the one-particle terms in that region, due to 
singularities in the rapidity variables. Hence we found that they can only contribute to the 
exponential factor. 

It would be very interesting to use these techniques to obtain the large-time expansion in differ- 
ent regions, or to solve numerically the linear integral equations in order to obtain correlation 
functions at all times. This would clarify the various approximate techniques that have been 
used up to now for studying finite-temperature Ising correlation functions, some of them being 
also applied to interacting integrable models. Our methods can be applied whenever there ex- 
ists an integrable PDE describing correlation functions (as is known to happen for twist fields 
in free fermionic models), and applying it to other models may further clarify the structure 
of finite-temperature correlation functions. Also, it could be instructive for the theory of in- 
tegrable PDEs to further investigate the connection between massive free models of QFT and 
integrable PDEs, and in particular, between the form factors and solutions to the GLM equa- 
tions. In particular, it would be interesting to understand further the origin of the complex 
region of the sinh-Gordon solution from the viewpoint of the PDE. For instance, it seems, al- 
though it is not entirely clear, that the information about how to go through the light-cone was 
already included in the initial scattering data, since the derivation of our solution did not make 
use of any asumption for crossing the light-cone. In a related direction, we wish to point out 
that, in a similar fashion to the zero-temperature case and the relation to Painleve equations, 
QFT ideas and results can give rise to non-trivial conjectures concerning solutions to integrable 
PDEs in the finite-temperature case. Finally, of course, our method cannot be directly applied 
to correlation functions in more complicated interacting integrable models, since there is gener- 
ically no PDE description. However, it may be worth looking for a generalisation of the GLM 
equations that give rise to form factor expansions at zero temperature in interacting integrable 
models; that may open the way to finite-temperature form factors in such models, to new ways 
of generating form factors, and to new representations for correlation functions. 
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A Derivation of the symmetric solution 



In this appendix, we show that certain objects in the QFT model satisfy the linear problem 
associated to the sinh-Gordon differential equation. Following [13] , let us consider the doubled 
theory: two copies (labelled by a and b) of the Majorana fermion, with anticommuting funda- 
mental fields and factorizing correlation functions. Symmetries in the doubled theory lead to 
trace identities for correlation functions incorporating the conserved charge. For example, the 
following integral is conserved: 



Y 1 



1 

2^ 

ni 

~2~ 



Tfl — 

4 J a{s)d s 4j b {s) + i—4j a {s)tp b {s] 



dx 



{A(pfB{0)+B{0fA(p)) 



d(3 
2^ 



where A(j3), B(j3) are the annihilation operators of particles with rapidity (3 for the copies a, b 
respectively (their hermitian conjugates are the corresponding creation operators) and ipai^b 
are the fundamental fermionic fields for copies a, b. The symbol s represents the space-time 
doublet (x, t), and the derivative d s = (d x — d t )/2 is the "holomorphic" derivative d that defined 
after (jSJ). Therefore, the following trace identity holds: 



[[Y 1 ,c a ( S y b (s)A\p')<T a (s')(j. b (s')]) 



T 



. 



(132) 



The twist fields a(s) are bosonic, whilst the /i(s) are fermionic. The commutator of the funda- 
mental field with twist fields may be derived from the definition of twist fields as highest weight 
states, see [13]. Those of interest to us are 



[Y 1 ,a a (s)a b (s)} 
[Yx,a a (s)ii b (s)) 



-d s n a (s)(i b (x) + fi a (s)d s ii b (s) 
-id s a a (s)fi b (x) + ia a (s)d s ^ b (s) 



(133) 



The trace identity defined above is therefore equivalent to 

< [- (d s fi a (s)) fi b (s) + /j, a (s)d s n b (s)\ A j (f3')a a (s')fi b (s')) T 

+ {a a (s)a b (s)A ] ((3') [i (d a > fJ, a (s')) cr b (s') - iji a {s')d sl a b {s')\ ) T = 



a a (s')n b (s')) T 
. (134) 



The correlation functions involving a product fields in both copies of the theory factorize. We 
define 

E(s + s', (3)F(s - s', (3) = (^s)A\p)a(s')) T 
E(s + s', (3)F(s - s', (3) ee (a(s)A\/3)fi(s')) T . 
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E(s,0) are the centre of mass plane waves e ip i 3X ^ 2 iE p 1 I 2 and the functions F are the finite 
temperature traces defined in section I3TT1 With these definitions, the above equation becomes 

d s (E(s + s')F{s - s'))G{s - s') - E{s + s')F{s - s')d s G(s - s') 

777 / 

+ °a(s)a b (s)— J e p B(J3y6(J3 - (3')a a (s')fi b (s')d(3 

+ id s ,{E{s + s')F(s - s'))G{s - s') - iE(s + s')F(s - s')d s ,G(s - s') = , 

where G, G are the correlation function of twist fields, as in the text. This expression may be 
simplified to (dropping the functional dependencies for clarity) 

d s (EF)G - EFd s G + GEF^- + id s >(EF)G - iEFd s ,G = . (135) 

The trace of the same conserved charge with a a (s)fi b (s)A^a a (s)a b (s') leads to a similar differ- 
ential equation 

ff 

~ ~ ~ 7YIP ~ ~ 

id s {EF)G - iEFd s G + GEF— d s ,(EF)G + EFd s ,G = . (136) 

Using: 

d s G= d -^G + d -^G 
2 2 

d s G= d -fG+ 9 -fG , (137) 
and also d s /G = —d s G, and d s E = d s /E = ime^E/A, and defining / = e~ x ^ 2 F, we find 

d s (f + if) = -i^(f - if) + d -f{f + if) 

d s (f - if) = -i^-U + if) ~ d -fU- if) ■ (138) 
A related set of equations are obtained from the conserved charge: 

Y-i = — / \^> a (s)d s ij; b (s) - i—ip b (s)il) a (s)j dx 



y f ^ (A(PYB(P) + B{ff? A{ft)) ^ 



where d s = (d x + d t )/2. A pair of relations like those in equation (11331) for Y-i, but with 
the replacements d s — > d s and i — >• — i on the right hand side, is used to obtain the following 
relations: 

Uf + if) = i r ^ 1 {f ~ if) ~ 8 -f{f + if) 

W ~ if) = i r ^-^U + if) + d ~fU- if) ■ (139) 

Using d x = d s + d s , a linear combination of equations (I138P and (I139P shows that the vector 
^ = (f — if, f+ if) satisfies the equation (d x — A x )^> = 0, where A x is the matrix associated 
to the zero-curvature formulation of the sinh-Gordon equation, as claimed in the text. It is 
referred to as the 'special' solution. 
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B Asymptotics for large and small A 



Here, we show how to extract the behaviour of a(A) and b(X, 0) as 9 — > ±00 by demonstrating 
the form of the Jost solution in these limits. This is one of the constraints on the form of a (A) 
in section 13.21 We find that the Jost solutions \P + and have the following A asymptotics for 
all x: 



e ip\x/2 I e -tpcr z /2 



e -ip\x/2 I e ~<pa z /2 



e ip\x/2 I e <p<r z /2 



-ip\x/2 I <pa z /2 



f- O 



|A| 



, |A| -> 00 



0(\X\) 



IAI ^0 



1 +0(|A|) , |A|-0 



(140) 
(141) 
(142) 
(143) 



The derivation of the first equation proceeds as follows. Start from the gauge transformed 
function 

\jr + = e i P\ x / 2 -^z/2-qj^ 
which satisfies the linear equation [d x + ip\(l — cr x )/2 — A x ]ty + = with 



if 2i(dt-d x )ip mX-^l-e 2 *?) 



4 V mX- x {l - e- 2t ?) -2i(d t - d x )y 



The boundary conditions on tp ensure that A x vanishes rapidly as x — > ±00. This differential 
equation may be rewritten as an integral equation: 

$l(x,\) = ( \ \ - £ Aye^-^-^A x {y,X)^{y,X). (144) 

Decomposing this function into the two eigenvectors of (1 — a x ), 

*l(x, A) = F(x, X) ( _\ ) + G(x, A) ( \ 



yields the coupled integral equations: 



F(x,X) = dye 



jp\{y-x) 



(d t - d x )ip 



G(y,X) 



im 



f — ( (1 - cosh(2^)) F(y, A) - sinh(2^)G(y, A) 



G(x, A) = 1 



dye 



im 
4A 



ip\(y-x) 



(d t - d x )<p 



F(y,X) 



[1 - cosh(2< / ?)) G(y, A) - smh(2<p)F(y, A) 



(145) 
(146) 
(147) 
(148) 
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Note that for f(x) decreasing sufficiently quickly as x — > oo 

POO 

/ dye lX ^f(y)=0(\- 1 ) , 

J X 

so that G(x, A) = 1 + 0(A _1 ) and F(x, A) = 0(A _1 ), as in equation ( 11401) . Derivations of the 
other three equations follow similarly. 

C Deriving the large-t expansion of the solution to the 
GLM equations 

In order to show that the asymptotic form of the solution to (11121) is, to some extent, unique, 
it is useful to derive it from a more general trial solution. Let us consider (I112p without the 
integral over F_iW\. The pole as x — > 2t of the function F^(x,t) makes solving the equation 
(11121) quite subtle: this pole may give a pole to y,t) as y — > x, which makes the integrals 
logarithmically divergent at the boundary z — x. Trying to take £(x,y,t) as a principal part 
with respect to x = y means that this divergency would be translated into some term log(e) - 
the integral would not be a true principal value integral, and £,(x,y,t) not a true distribution. 
Due to the term with 6(2t — x — y), there is also a possible problem at y = 2t — x in y, t), 
where there could be a finite jump (in addition to a pole). To make everything clear, we need 
to regularise the GLM equations. One could take from the beginning fl98l) with finite e for the 
integrals (11011) . or one could keep the number e of the principal part prescription for F^(x,t) 
finite in the integral equations. For our purposes, it will be clearer to modify the integral over 
z in the GLM equation (11121) by requiring z — x > e and \z — (2t — x)\ > e for some fixed 
positive number e. Then, we may see the solution as a function of e, and the integral equations 
are true principal value integrals in the integration regions that are not excluded. We seek a 
solution that gives finite values for the quantities associated to the PDE, ^± and ip, as e — ■> 0; 
we will find a solution that is itself finite as e — ► 0. 

It will be convenient to replace the variables x, y with A = — x + t and s = y — t. Our large 
t region corresponds to A,m _1 ,T _1 <^ t. Also, from the discussion above, the solution £ may 
separate into two regions (between which there could be a discontinuity): 

• "spacelike", £ S (A, s,t): s > A 

• "timelike", f*(A, s,t): -A < s < A . 

Discarding for a moment the contribution from the integral over F_iW\ (since it will be shown 
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to contribute subleading terms), in terms of these variables, the GLM equation becomes 
- J^£'(A, s , t) = F p (2t + s - A, t) ( | ) + F p (2t + s + A, t) ( ~? ) 



/A — £ /"OO 
F p (2t + s + s', t)f'(A, s', *)ds' + / F p (2t + s + s', t)f (A, s', i)da' 
-A+£ </A+£ 

+ -e < (A,-s,t) (149) 
m 



and 



(A, s, = F P (2t + s - A, t) ( \ ) + F p (2t + s + A, i) ( . l ^ 

/A — £ /"OO 
F p (2t + s + S , ,t)£'(A,s / ,t)ds , + / F p (2* + s + sV)£ s (A,s',f)ds' . (150) 
-A+£ J A+e 



A linear combination of these equations and the same equations with s — > —s yields the coupled 
equations for £ s and £*: 



^ (A ' M) = ~ f " F p (2t + s + s>, t)e'(A, S ', t)d S ' - ^ r F p (2t + S + a ', t)r(A, s', t)ds' 
m 4 J_ A+e 4 J A+£ 

■ pA-e ■ roo 

+ - / F p (2t -s + s', t)£ 4 (A, s ', t)ds' + - / F p (2t -s + s', t)£ s (A, s', t)ds' 
+ ^F p (2t + s + A,t) ( M + ^F p (2t - a + A,t) ( J ) 

- ^F p (2t + 8 - A, t) ( ) ) + ^F p (2t - S - A, t) ( J ) , (151) 
and 

P(A,M) = °£ /" A ^ Fo p (2t + s + S ',t)^(A, S ',t)d S '-^ r F p (2t + s + s',t)C(A,s',t)ds' 

m l J-A+e 1 J A+e 

+ V p (2t + S + A,t)( J ) -?p p (2t + s -A,t)( | ) . (152) 

We are looking for the leading terms of the solution to these equations. More precisely, we will 
concentrate on terms which, in the limit t — > oo with s, A fixed, behave like e ~ l W^-\/l s ± A l and 
e -my/2i^y\s±A\ ^ ^ mes p 0wers Q f f_ For instance, choosing s = 0, this gives a precise large-time 
frequency of oscillation or exponential decay, and products of such terms give higher frequencies 
or stronger exponential decays. Terms with higher exponential decay or higher frequency will 
be deemed sub-leading. Hence, we will neglect all terms with v ^ in (11171) (and only the first 
series, with v — 0, remains). It is these leading terms that will give the leading terms in our 
region of interest for the function tp itself. With that understanding, for ease of notation, let 
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us define the following four functions: 

q(A,s,t) =F[(2t + s + A,t) 
u(A, s, t) = F p (2t - s + A, t) 
v(A, s, t) = F p (2t + s - A, t) 
iu(A, s, t) ee F p (2t - s - A, t) . 



(153) 



These are the inhomogeneous parts of the GLM equations for (equation 11511) and £ s (equa- 
tion [152]). Let us consider the trial solutions 



= (ag + /3m + 'jv + 5u>) 



f = ( a 'q + (3'u + yu + 8'w) 



(154) 



where the coefficients are two by two matrices in R + a z M, and insert them into ( 1151ft and 
( 11521) . It turns out that the leading terms (in the sense above) resulting from the integrals 
are of similar functional form, namely t- and e-dependent linear combinations of q,u,v,w and 
of their s-derivatives to arbitrary high orders. Some details of the calculations to obtain the 
leading terms are done below. In fact, the leading power of t° in these linear combinations does 
not contain any derivative, so that to leading order these trial solutions give an inhomogeneous 
linear matrix equation for the coefficients. In matrix form, this equation is, to 0(e), 



M 



( 


1 

2 




ia z 
2 


<J Z K 
4 








a z K 
4 


o z K 
4 


o\ 




ia z 


iK 


1 


, iK 
4 








iK 


iK 







2 


4 


2 


4 


1 




o z K 
4 







1 

2 


ia z i (T Z K 
2 ' 4 
















iK 
4 







irr z 
2 


1 iK 

2 4 



















-io z - ^ 


-ia z 


a z K 
2 





1 


cr z K 
2 


cr z K 
2 
























1 










2 










a z k 

2 








1 





{ 

























V 



M 



(a\ 

7 
5 



m 



( i/4 \ 

^/4 
-^/4 
i/4 
i/2 


-a z /2 




a 
0' 

i 

47 + 41n(mv / it/2), with 7 Euler's constant. The linear 
combinations to all orders in negative powers of t are formally reproduced by this matrix if 
K and K are, instead, power series in t™ 1 whose coefficients are differential operators in s of 
order equal to the power of t~ l (and then, K + K 7^ — 2i). The matrix above is invertible for 
arbitrary K and K and has the miraculously simple solution: 



in which K + K = —2i and ttK 



ma. 



mcr z 



(155) 
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Since this is independent of K and K, it is correct to 0(t 00 ). This gives the solution written 
in the text, ffTTBD . 

In order to obtain the matrix above, we use the expansion (TTTTfl for the function F p (x,t) (the 
part with v = only). The idea is that if the large-t limit (with A fix) can be taken in the 
integrands of (115ip . (I152p . then the result is sub-leading. The only problem in taking this 
limit is at the points in the integration region where the coefficient of t in the argument of the 
Bessel functions in (I117P vanishes. This happens only for 2x = t in F^(x,t) (that is, at its 
singular point - this is why the regularisation e is important). There are two types of points 
where this happens. First, the kernels F^(2t ± s + s') of the integral equations (11511) . (11521) 
give a principal-value integral about the point ±s + s' = 0. The s'-contour can completed and 
deformed towards the positive imaginary direction (this direction is chosen due to the square- 
root structure of the argument of the Bessel functions), and the half-residue substracted. The 
resulting half-residues give contributions to the leading terms of interest. 

The second type of points are those coming not from the kernels, but from the functions q, u, v, w 
in the trial solution. For the functions q and w, they are at the boundary s' = —A + e; for u 
and v, they are at A ± e. Let us give an example with the function q. We have, for arbitrary 

g(s'), 

[ g(s')q(A,s',t)ds' = [ g(s')F p (2t + s' + A)ds' 

J-A+e J-A+e 

~ l 9 ^'-^ (i&rh)) K i-^ m VK2t + s'/2)(-is'/2))ds' 

With the change of variable in the last line, the point s' = —A has been "resolved", and the 
large-t expansion can be taken in the integrand. This gives 

m- l {Kg){-A) 

where K is an expansion in inverse powers of t with coefficients that are, generically, differential 
operators. More precisely, 

K = 4 7 + 4 ln(m v / it/2) + —(^- + . . . 

mt \ 7r Tim J 

where s is the independent variable of the function on which it applies. Similarly, 

/ g(s')w(A,s',t)ds' ~m-\Kg){-A) 

J-A+e 
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where K is a different operator. Other integrals give 

•A-e 



g(s')u(A,s',t)ds' ~ m~\{-2i-K)g){A) 

g(s')v(A,s',t)ds' ~ m~\(-2i-K)g)(A) 

/ g(s')u(A,s',t)ds' ~ m" 1 (^)(A) 

JA+e 

[ g(s')v(A,s',t)ds' ~ m" 1 (A^)(A). 

./A+e 

Using these, one gets the other terms giving the 8 by 8 matrix above. 
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